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Abstract. Since 1997 tlie neutrino telescope 
AMANDA at the geographic South Pole has been 
monitoring our Galaxy for neutrino bursts from su- 
pernovae. Triggers were introduced in 2004 to submit 
burst candidates to the Supernova Early Warning 
System SNEWS. From 2007 the burst search was 
extended to the much larger IceCube telescope, which 
now supersedes AMANDA. By exploiting the low 
photomultiplier noise in the antarctic ice (on average 
280 Hz for IceCube), neutrino bursts from nearby 
supernovae can be identified by the induced collective 
rise in the pulse rates. Although only a counting 
experiment, IceCube will provide the world's most 
precise measurement of the time profile of a neutrino 
burst near the galactic center. The sensitivity to 
neutrino properties such as the 6*13 mixing angle and 
the neutrino hierarchy are discussed as well as the 
possibility to detect the deleptonization burst. 

Keywords: supernova neutrino IceCube 

I. Introduction 

Up to now, the only detected extra-terrestrial sources 
of neutrinos are the sun and supernova SN1987A. To 
extend the search to TeV energies and above, neutrino 
telescopes such as AMANDA and IceCube 1 1 1 have been 
built. It turns out that the noise rates of the light sensors 
(OMs) in the antarctic ice are very low 700 Hz 
for AMANDA, ~ 280 Hz for IceCube) opening up 
the possibility to detect MeV electron anti-neutrinos 
from close supernovae by an increase in the collective 
rate of all light sensors. The possibility to monitor the 
galaxy for supernova with neutrino telescopes such as 
AMANDA has first been proposed in 121 and a first 
search has been performed using data from the years 
1997 and 1998 (1. 

The version of the the supernova data acquisition 
(SNDAq) covered in this paper has been introduced for 
AMANDA in the beginning of the year 2000 and was 
extended to IceCube in 2007. The AMANDA SNDAq 
has been switched off in February 2009. We will in- 
vestigate data recorded by both telescopes concentrating 
on the 9 years of AMANDA measurements and make 
predictions for the expected sensitivity of IceCube. 

II. Detectors and Data Acquisition 

In AMANDA, the pulses of the 677 OMs are collected 
in a VME/Linux based data acquisition system which 
operates independently of the main data acquisition 
aimed at high energy neutrinos. It counts pulses from 
every connected optical module in a 20 bit counter in 
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Fig. 1: Rate distribution of a typical IceCube module 
The recorded rate distribution (filled area) encompasses about 44 
days and has been fitted with a Gaussian (solid line, x^/^doi = 
56.5) defined by fi = 271 Hz and cr = 21 Hz. However, a 
lognormal function (dotted line, x^/^doi = 0.66) with geomet- 
ric mean figco =6.41n(Hz), geometric standard deviation o-gco = 
0.031n(Hz) and a shift of xq = 377 Hz is much better at describing 
the data. 



fixed 10 ms time intervals that are synchronized by a 
GPS-clock. 

In IceCube, PMT rates are recorded in a 1.6384 ms 
binning by scalers on each optical module. The infor- 
mation is locally buffered and read out by the IceCube 
data acquisition system. It then transfers this data to the 
SNDAq, which synchronizes and regroups the informa- 
tion in 2 ms bins. 

The software used for data acquisition and analysis 
is essentially the same for AMANDA and IceCube. The 
data is rebinned in 500 ms intervals and subjected to an 
online analysis described later. In case of a significant 
rate increase ("supernova trigger"), an alarm is sent to 
the Supernova Early Warning System (SNEWS, |4|) via 
the Iridium satellite network and the data is saved in a 
fine time binning (10 ms for AMANDA and 2 ms for 
IceCube). 

HI. Sensor Rates 

In 500 ms time binning, the pulse distribution of 
the average AMANDA or IceCube OM conforms only 
approximately to a Gaussian. It can more accurately be 
described by a lognormal distribution (see figure [TJ. 

The pulse distributions exhibit Poissonian and cor- 
related afterpulse components contributing with similar 
strengths. The correlated component is anticorrelated 
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with temperature and arises form Cherenkov light caused 
by decays and glass luminescence from radioactive 
decay chains. A cut at 250 /iS on the time difference 
of consecutive pulses effectively suppresses afterpulse 
trains, improves the significance of a simulated super- 
nova at 7.5 kpc by approximately 20% and makes the 
pulse distribution more Poissonian in nature. 

IV. Effective Volumes 

Supernovae radiate all neutrino flavors, but due to the 
relatively large inverse beta decay cross section, the main 
signal in IceCube is induced by electron anti-neutrinos 
(see Q). The rate R per OM can be approximated 
weighing the energy dependent anti-electron neutrino 
flux at the earth (derived from the neutrino 

luminosity and spectra found in f6^|) with the effective 
area for anti-electron neutrino detection Aos{Ep^) and 
integrating over the whole energy range: 

J^dE^^ ^{E^JA^siE^J , with 
Aeff (E^^ ) = n j;^dE,+ ^ {E,^ , E,+ ) V,s,e^ {E,+ ) • 

dE'+ (-Ep^ , i?e+) is the inverse beta decay cross section, 
n the density of protons in the ice and Ygg e+ (£^e+ ) the 
effective volume for positron detection of a single OM. 

y^ff g+ (i?g+ ) can be calculated by multiplying the 
number of Cherenkov photons produced with the effec- 
tive volume for photon detection Vcff,'y„h- 

By tracking Cherenkov photons in the antarctic ice 
around the IceCube light sensors [71 and simulating 
the module response one obtains V^ff.-y^j, — 0.104m'^ 
for the most common AMANDA sensors and Vqb,^ = 
0.182 for the IceCube sensors. With a GEANT-4 
simulation, the amount of photons produced by a 
positron of an energy E^+ can is estimated to be 
iVT"^!' = 270£;e+/McV. Consequently, the effective 
volumes for positrons as a function of their energies 
are ^^^^^^ - 19.5 E,^^ and V^^f^^'^'^ = 

34.2i?g+ j^y. Uncertainties in the effective volumes 
derive directly from the uncertainties of the ice models 
(-- 5%) and in the OM sensitivities (-- 10%). 

In this paper we assume a supernova neutrino pro- 
duction according to the Lawrence-Live nno re model 
as it is the only one that provides spectra for up to 
15 s. It gives a mean electron anti-neutrino energy of 
15 MeV corresponding to an average positron energy of 
(13.4 ± 0.5) MeV. 

V. Analysis Procedure 

A simple investigation of the rate sums would be 
very susceptible to fluctuations due to variations in the 
detector response or external influences such as the 
seasonal variation of muon rates. Medium and long term 
fluctuations are tracked by estimating the average count 
rate by a sliding time window. The rate deviation for 
a collective homogeneous neutrino induced ice illumi- 
nation is calculated by a likelihood technique. In time 
bins of 0.5 s and longer, the pulse distributions can be 



approximated by Gaussian distributions. For a collective 
rate increase A^, the expectation value of the average 
mean rate /i^ of a light sensor i with relative sensitivity 
ei increases to /i^ + e^A/i. The mean value /i^ and 
its standard deviation cr,; are averaged over a sliding 
window of 10 min, excluding 15 s before and after the 
0.5 s time frame studied. By taking the product of 
the corresponding Gaussian distributions the following 
likelihood for a rate deviation Afi is obtained: 



C 



NoM 

n 

i=l 



1 



'271 a. 



exp - 



2a 




Minimization of — ln£ leads to: 
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The data is analyzed in the three time binnings 0.5 s, 4 s 
and 10 s for the following reasons: First, the finest time 
binning accessible to the online analysis is 0.5 s. Second, 
as argued in |9|, the neutrino measurements of SN1987A 
are roughly compatible with an exponential decay of 
T = 3s. The optimal time frame for the detection of 
a signal with such a signature is w 3.8 s. Last, 10 s are 
the approximate time frame where most of the neutrinos 
from SN1987A fell. 

To ensure data quality, the optical modules are sub- 
jected to careful quality checks and cleaning. Those 
modules with rates outside of a predefined range, a high 
dispersion w.r.t. the Poissonian expectation or a large 
skewness are disqualified in real time. 

Since SNEWS requests one alarm per 10 days, the 
supernova trigger is set to 6.3(7. 

To ensure that the observed rate deviation is homo- 
geneous and isotropic, the following discriminant is 
examined: 
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We demand the data to conform to a x'^ -confidence 
level of 99.9%. However, it was found that the 
cannot clearly distinguish between isotropic rate changes 
and fluctuations of a significant number of OMs: If 
e.g. 20% of the OMs record rate increases of about 
1 cr (« 20 pulses/0.5 s), the significance for isotropic 
illumination can rise above 6 a without being rejected 
by the condition. Still, no method was found that 
performed better than x^- 

VI. External Perturbations 

Figure |2] shows the distribution of the significance ^ — 
for AMANDA. From the central limit theorem, one 
would expect a Gaussian distribution with a width of 
a ~ 1. This expectation is supported by a background 
simulation using lognormal representations of individual 
light sensor pulse distributions. However, one finds that 
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Fig. 2: Significances of AMANDA in 0.5 s time frames 
The filled area shows the significance of the data taken in 2002 with 
a spread of a = 1.13. The solid line denotes the initial background 
simulation (cr « 1) and the dotted line is a toy monte carlo taking 
into account the fluctuations due to muon rates (cr « 1.11). 



the observable is spread wider than expected and exhibits 
a minor shoulder at high significances. 

As this observable is central to the identification of 
supernovae, its detailed understanding is imperative and 
a close investigation of these effects is necessary. 

Faults in the software were ruled out by simulating 
data at the most basic level, before entering the SNDAq. 
Hardware faults are improbable, as the broadening of 
the rates is seen both for AMANDA (cr = 1.13) and 
IceCube (a = 1.27) and the detectors use different and 
independent power and readout electronics. 

One can then check for external sources of rate 
changes as tracked by magnetometers, rio- and pho- 
tometers as well as seismometers at Pole. All have 
been synchronized with the rate measurements of the 
OMs. Only magnetic field variations show a slight, 
albeit insignificant, influence on the rate deviation of 
—4- 10~^ S|. for AMANDA. Due to a /i metal wire mesh 
shielding, the influence on IceCube sensors is smaller by 
a factor ~ 30. 

It turns out that the main reason for the broadening 
are fluctuations of the atmospheric muon rates. During 
0.5 s, AMANDA detects between w 3.0- 10^ (June) and 
« 3.4-10'^ (December) sensor hits due to muons. Adding 
hits from atmospheric muons broadens and distorts the 
distribution derived from noise rates. A simulation taking 
into account the muon triggered PMT hits increases the 
width of the significance distribution to « 1.11 (see 
figure |2]i. Although the investigations are still ongoing 
for IceCube, we expect the larger a can be ascribed to 
its lower noise rate leading to a higher resolution and 
thereby a stronger sensitivity to perturbations. 

A Fourier transformation was performed on the data 
stream. No evidence for periodically recurring events 
was found, but occasional correlations between subse- 
quent 0.5 s time frames both in the summed noise rate 



and the rate deviation could be identified. We determined 
the mean significance before and after a rate increase 
exceeding a predefined level. One observes a symmetric 
correlation in the data mainly between ±10 s. The origin 
of the effect, which is present both in AMANDA and 
IceCube data, is under investigation. 

VII. Expected Signal and Visibility Range 

A preliminary analysis of AMANDA data from 2000 
to 2003 yielded a detection range of R^^ = 14.5 kpc 
at the optimal binning of 4 s at an efficiency of 90%, 
encompassing 81% of the stars of our Galaxy. As signals 
with exponentially decreasing luminosity (at r = 3 s) 
have been used as underlying models, the two other 
binnings were less efficient with i?o.5s ~ 10 kpc (56%) 
and i?ios = 13 kpc (75%), respectively. 

The rates of the IceCube sensors are stable and 
uniform across the detector. Scaling the observed rates 
to the full 80 string detector, a summed rate of 
(i?iccCubo) = (1.3-10^±1.8-102)Hz is expected. While 
a single DOM would only see an average rate increase of 
13 Hz or 0.65 a, the signal in the whole IceCube detector 
would be 6.1 • 10** Hz or 34 a for a supernova at 7.5 kpc 
distance. Using this simple counting method, IceCube 
would see a supernova in the Magellanic Cloud with 
5 cr significance. 

Figure [3] shows the expected signal in IceCube for 
a supernova at 7.5 kPc conforming to the Lawrence- 
Livermore model with w 10^ registered neutrinos in 
15 s and a statistical accuracy of 0.1% in the first 2 s. 
Assuming 2 • 10^ events in Super-Kamiokande (scaled 
from 1 10 1), one arrives at an accuracy of 1% in the same 
time frame. While IceCube can neither determine the 
directions nor the energies of the neutrinos, it will pro- 
vide worlds best statistical accuracy to follow details of 
the neutrino light curve. Its performance in this respect 
will be in the same order as proposed megaton proton 
decay and supernova search experiments. As the signal 
is seen on top of background noise, the measurement 
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Fig. 3: Expected signal of a supernova at 7.5 kpc in IceCube. 
The solid line denotes the signal expectation and the bars a randomized 
simulation. 
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Fig. 4: Supernova detection ranges using 0.5 s time frames 
Significances which supernovae conforming to the Lawrence- 
Livermore model would cause as function of their distance. IceCube 
performance is described by the solid line, AMANDA by the dotted 
Hne. Both are given for the 0.5 s binning. 

accuracy drops rapidly with distance. Figure |4] shows the 
significance at which IceCube and AMANDA would be 
able to detect supernovae. 

VIII. Deleptonization Peak and Neutrino 
Oscillations 

As mentioned before, the signal seen by AMANDA 
and IceCube is mostly dominated by electron anti- 
neutrinos with a small contribution by electron neutrino 
scattering, thereby leading to a sensitivity which is 
strongly dependent on the neutrino flavor and is thus 
sensitive to neutrino oscillations. 

As the neutrinos pass varying levels of density within 
the supernova, the flux of electron and electron anti- 
neutrinos (1> released is different from the initial flux $° 
produced during the collapse: 

= + ■ 

The survival probability p/p for the v^/ve^ varies with 
the mass hierarchy and the mixing angle 6*13 ifTTl : 
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At the onset of the supernova neutrino burst during the 
prompt shock, a ^ 10 ms long burst of electron neutrinos 
gets emitted when the neutron star forms. As the shape 
and rate of this burst is roughly independent of the 
properties of the progenitor stars, it is considered as 
a standard candle, allowing one to determine neutrino 
properties without knowing details of the core collapse. 

The first 0.7 s of a supernova signal were modeled af- 
ter |[T2l . Figure |5] shows the expectation for a supernova 
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Fig. 5: Neutrino signals of a supernova at 7.5 kpc distance modified 
by oscillations as seen in IceCube. 

The line shows the expectation without neutrino oscillations, the 
squares show the simulated signal for the normal mass hierarchy 
with sin^ S13 > 10"'' (sin^ 613 < 10~^ lies nearly on the same 
line) and the triangles show the signal for inverted mass hierarchy at 
sin2 6»i3 > 10"^. 

at a distance of 7.5 kpc in the 2 ms binning of IceCube. 
Due to statistics and the rising background from the 
starting electron anti-neutrino signal, the identification 
of the neutronization burst is unlikely at this distance. 
However, with a supernova at 7.5 kpc one might be able 
to draw conclusions for the mass hierarchy, depending 
on the reliability of the models. 

IX. Conclusions and Outlook 

Wifli 59 strings and 3540 OMs installed, IceCube 
has reached 86% of its final sensitivity for supernova 
detection. It now supersedes AMANDA in the SNEWS 
network. With the low energy extension DeepCore, 
IceCube gains 360 additional OMs with a ~ 30% higher 
quantum efficiency (rates ^ 380 Hz). These modules 
have not been considered in this paper. 
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